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THE  k  SHORTEST  ROUTES  AND  THE  k  SHORTEST  CHAINS  IN  A  GRAPH 


by 

Michel  Sakarovitch 


1.  Introduction 

In  a  given  graph  G  =  (N,A)  where  N  represents  the  set  of  nodes  and  A  the 

set  of  directed  arcs,  a  length  a(x,y)  is  associated  to  every  arc  (x,v)  e  A  .  A 

route  from  x  to  x  is  a  sequence  of  nodes  and  arcs 
on 


P  =  [xq, (xq,x1) ,x1§ (x1#x2) 


x  . , (x  ,  ,x  ) ,x  ] 
n-1  n-1  n  *  n 


such  that 


x^  e  N 

(Vi.*i)  e  A 


i  0 , « « « «  ,n 
i  e  1 . n 


A  route  P  will  also  be  considered  as  a  partial  subgraph  (see  [2]  )  of  G 

containing  nodes  x  ,..,x  and  arcs  (x  .x_)  .....  (x  ,  ,x  )  .  The  length  of  a 

o  n  o  1  n-1  n  ° 

route  P  is: 

n 

£(P)  -  l  a(x  ,x  ) 
i=l  1-1 


A  chain  is  a  route  in  which  no  node  is  repeated.  A  cycle  is  a  route  in  which 


no  node  is  repeated  except  that  xq  =  x^  .  0  and  D  will  be  two  particular  nodes. 
P  (P1)  will  denote  the  shortest  route  (chain)  from  0  to  D  and: 

q  q 


G 

q 


<Nn»An> 

q  q 


q 

u  P 
J-i 


j 


will  be  the  partial  subgraph  of  G 


x  e  N 

iff 

q 

(x,y)  e  Aq 

iff 

defined  by: 


x  e  Pj  for  some  j  £  q 
(x,y)  e  P.  for  some  j  <_  q 

J  • 
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In  a  similar  way  G  * ,  A',  A'  and  N'  are  defined,  k  is  the  number  of  routes 

q  q  q  q  q 

in  G  which  we  need  to  generate  as  P^..  .  Pq„  i.e,  k  is  the  min  of 

q  q+i  q+kq  q 

(k  -  q)  and  (max  number  of  routes  in  G  -  q)  . 

Our  central  problem  will  be  the  construction  of  the  k  shortest  routes  (k 
shortest  chains)  from  0  to  D  .  In  this  form,  the  problem  is  slightly  more 
general  than  the  one  dealt  with  in  [1],  [6]  and  [3]  and  reviewed  in  [7]  in  two 
respects: 

-  Negative  lengths  are  allowed,  provided  that  no:  "negative  cycle"  (a  cycle 
along  which  £  a(x,y)  <  0)  exist. 

-  Attention  is  not  restricted  to  chains. 

Moreover,  this  paper  is  based  on  the  following  simple  idea: 

At  any  stage  of  the  algorithm,  suppose  we  know  the  q  (q  <  k)  shortest  routes 
from  0  to  D  and  p|j+1  Pq+k  which  are  the  (q+l)St  .  (q+k^)th 

shortest  routes  in  ordered  according  to  increasing  lengths.  is  a  shortest 

■ ■ ■  — 

route  from  0  to  D  which  has  at  least  one  arc  in  A  .  The  (q+1)  shortest 

q 

route  is  obtained  by: 


p’+i  if  5 

Pq+1  “  l(Pq+lJ  <  *<P’+1> 


G^  and  the  list  P^j+^  *s  updated  and  the  process  begins  again. 

The  algorithm  is  initiated  by  finding  the  shortest  routes  from  0  to  every 
node  using  an  algorithm  such  as  the  one  in  [A]  if  all  lengths  are  nonnegative  and 
the  one.  in  [5]  if  there  are  negative  lengths. 


In  the  next  section  an  algorithm  to  construct  the  k  shortest  routes  from  0 


to  D  is  described  and  it  is  justified  in  Section  3.  Since  the  problem  of  finding 
the  k  shortest  chains  is  of  interest  in  many  applications,  a  variant  of  the 
algorithm  presented  in  Section  2  which  exhibits  the  k  shortest  chains  is  described 
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in  Section  4.  Attention  is  drawn  to  the  fact  that  minor  modifications  of  these  two 
algorithms  allow  to  find  the  routes  (chains)  shorter  than  a  given  length.  If 
capacities  were  to  be  considered,  an  immediate  application  of  this  last  problem 
could  be:  what  quantity  can  be  shipped  from  0  to  D  if  we  demand  that  the 
maximum  time  of  travel  be  less  than  or  equal  to  a  certain  given  number. 


2.  Algorithm  For  Finding  The  k  Shortest  Routes 

Step  0.  -  Find  all  the  shortest  ruutes  from  0  to  D  :  P  ,  ...,P 
If  q  ^  k  ,  terminate. 

Determine  for  all  x  e  N  : 

tt(x)  =  shortest  distance  from  0  to  x  . 
n(x)  =  shortest  distance  from  x  to  D  . 

Note  that  Pq+^  does  not  exist  (the  length  of  a  route  which  does 
not  exist  is  taken  as  infinite) . 
k  *=  0  . 

q 

Step  1.  -  Let: 


(1) 


Vx)  = 


Min 


(y,x)  e  A 


tt(x) 


for 

x  e 

N 

( 

for 

x  i 

N 

(2) 


E  **  {x  e  N  |  C  (x)  +  n(x)  =  Min  [c  (y)  +  n(y)]  =  n  ) 
q  yeN  q  1 


If  ir  =  00  and  k  +  q  <  k  :  terminate,  the  problem  is  infeasible. 

q  q 

m  by : 


*<Pq+m>  i  \  <  « W 


Determine 
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Let: 


G  —  G 
q+m  q 

pq+™  =  pq 

q+m+i  q+m+i 
q  -  q  +m 

j  =  0 


i  =  1, .  .  ,m 


i  =  1, , (k^-m) 


For  each  x  c  E  generate  all  the  routes  from  0  to 

x  ,  the  last  arc  of  which  belongs  to  ,  and,  the  length  of  which 

is  Cq(x)  >  this  is  easily  done  by  backtracking  in  the  definition  of 
C^(x)  and  n(x)  .  Similarly  generate  all  the  chains  from 

x  to  D  of  length  n(x)  . 

By  combining  a  route  R(ox)  with  a  chain  C(xD)  and  repeating 

the  operation  for  all  x  e  E  ,  all  the  routes  from  0  to  D  of 

length  rr  can  be  generated.  Lets  call  them  Pq  , ...  Pq,  .  We 

q  q+1  *  q+n 

need  only  to  generate  at  most  k  -  q  of  them. 


If  n  >  k  -  q  let  P 


«  Pq 


q+j  q+j 

solution  has  been  obtained. 


j  =  l,....,k  -  q  ;  terminate, 


Go  to  Step  2. 

Step  2.  -  Let  j  =  j  +  1 

If  j  >  n  go  to  Step  1, 


>  = 

q+i  q+j 


G  =  G  U  Pq 
q+1  q  q+j 
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If  the  set  of  arcs  in  G  ,  -  G  is  empty  let: 

q+i  q 


I  <  i  <  k 


pq+1  =  pq 

q+l+i  q+i 
q  =  q  +  1 

Go  to  Step  2. 


If  the  set  of  arcs  in  G^+^  -  is  not  empty,  they  form  a  simple 

route  C.  from  0.  to  D,  (0.,  D.  e  G  )  (of  Lemma  5). 

j  J  J  3  J  q 

Let  R  =  (R  ,..,R  )  ,  S  =  (S . Sg)  ,  T  =  (T  , . . ,T  )  be 

respectively  a  set  of  routes  from  0  to  0_.  ,  from  D_.  to  CL 
and  from  D_.  to  D  ,  obtained  by  considering  parts  of: 

-  Those  routes  from  0  to  D  which  contain  0  ; 

J 

-  Those  routes  from  0  to  D  which  contain  D.  and  then 

J  - 

°j  5 

-  Those  routes  from  0  to  D  which  contain  , 

Generate  new  routes  in  G^  ,  according  to  increasing  lengths, 

by  taking  the  union  of  a  part  of  a  route  in  R  ,  from  0  to  0^  , 
route  Cj  and  a  part  of  a  route  in  T  from  to  D  .  Insert 

those  newly  created  routes  in  the  sequence  so  that  the  new 

sequence  is  still  ordered  according  to  increasing  lengths.  Generate 
new  routes  until  exhaustion  of  possibilities  or  until  the 


Pq-k-q-n  +  j-1 


first  elements  of  this  new  sequence  have  been  created. 

If  S  =  0  go  to  Step  4. 

Let  u  =  1  .  We  will  call  route  of  order  u  a  route  which  contains 
(u+1)  times  route  (and  u  times  an  element  of  S  ).  For 

instance  a  route  of  order  1  is  the  union  of  a  route  in  R  ,  route 
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C.  ,  a  route  in  S  ,  route  C.  ,  a  route  in  T  . 
J  3 

Go  to  Step  3. 


Step  3.  -  Generate  routes  of  order  u  according  to  increasing  lengths  and 
their  elements  in  the  sequence  just  generated  in  such  a  position 
that  the  sequence  is  still  ordered  according  to  increasing  lengths 
either  until  exhaustion  of  possibilities  or  until  the  first 

elements  of  the  sequence  have  been  generated. 

If  the  first  route  of  order  u  is  among  the  p^  first  elements  let 
u  -  u  +  1  and  go  to  Step  3. 

n 

If  not,  go  to  Step  A. 

Step  4.  -  Let  q  =  q  +  1 

Denote  by  dements  of  the  sequence  just  generated  and  by 

k  their  number. 

q 

Go  to  Step  2. 

3.  Justification  Of  The  Algorithm 
Let 

A  “  { (x,y)  e  A  |  7r(y)  -  tt(x)  =  a(x,y)} 

G  =  (N,A) 

Lemma  1.  -  A  necessary  and  sufficient  condition  for  a  cycle  to  be  of  0  length 
(0-cycle)  is  that  it  be  contained  in  G  . 


Proof : 


(i)  If  a  cycle  [xo>x^, . . . »xn_^>xn  =  XQ]  belongs  to  G  ,  it  is  a  0  cycle. 


tt(x  ) 

i J 
I 

"(x  ) 
o 


-  k(x  )  =  a(>:  , x  ) 


c’  1 


Tf(x  J  =  a(x  > x  ) 
n-i  n-1  o 


=  l  n  a^xi_i>xi) 


0 


(ii)  Let  C  be  a  0  cycle.  For  all  the  arcs  (x,y)  e  C  we  have  by 


definition  of  if  : 


n(y)  -  ir(x)  <_  a(x,y) 

If  for  one  arc  we  had  strict  inequality  by  adding  these  relations  for 
all  arcs  of  C  ,  we  would  obtain  0  <  £(C)  which  contradicts  the 
assumption. 

q.e.d 

Lemma  2 .  -  All  routes  in  G  from  0  to  a  given  node  z  have  the  same  length: 
tt(z)  .  Any  route  not  completely  included  in  G  has  a  larger  length. 


Proof: 

Obvious,  since  along  any  route  of  G 

if(y)  -  "(x)  =  a(x,y) 

and  for  a  route  not  completely  included  in  G  ,  we  would  have  on  at  least  one  arc: 
Ti(y)  -  ti(x)  <  a(x,y)  . 

q.e.d 


Lemma  3.  -  A  necessary  and  sufficient  condition  for  n  to  be  finite  at  the  end  of 
Step  1  is  that  P^,.  do  not  contain  a  0-cycle. 

q+j  ' 


Proof : 


Suppose  P^+_j  contains  a  0-cycle,  then  the  routes  obtained  by  the  union  of 
P^+j  and  i  times  this  0-cycle  are  of  same  length  for  any  positive  integer 

i  . 


If  there  is  an  infinite  number  of  routes  of  length  ir^  ,  there  is  an  infinite 
subset  of  these  routes  the  elements  of  which  differ  from  one  another  by  the 
addition  of  a  cycle  a  certain  number  of  times  since  the  number  of  chains  from  0 


to  D  is  finite.  The  length  of  this  cycle  is  thus  0  . 


q  .e.d 


Lemma  4.  -  it  ,  as  defined  in 
-  q 

completely  included  in  G  . 

q 


(2)  ,  is  the  shortest  length  of  all  routes  not 


Proof: 

* 

It  is  sufficient  to  show  that  for  each  route  P  which  contains  at  least  an 

—  'k  m " 

arc  (x,y)  in  A  ,  £(P  )  >  it  . 

q  ’  —  q 

From  the  definition  of  ir  and  n  we  have: 

(3)  £(P  )  >_  tt(x)  +  a(x,y)  +  n(y) 


But  from  (1)  : 


tt(x)  +  a(x,y)  >_  ^^(y) 

and  the  property  of  shortest  distances  that: 

Tf(y)  <.  ti(x)  +  a(x,y) 


Then,  from  (2)  : 


4q(y)  +  n(y)  1  ^q 

which  together  with  (3)  implies: 

✓ 

£(P*)  >  ii 

~  q 


q.e.d 


Lemma  5.  -  In  Step  2  the  set  of  arcs  in 

simple  route  C.  from  0.  to  D.  (0., 
j  j  3  3 

or  a  simple  cycle  in  which  case  D.  =  0. 

J  J 


G  , ,  -  G 
q+1  q 

D.  c  G  )  . 

J  q 


is  either  empty  or  forms  a 
This  route  is  either  a  chain 
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Proof: 

If  the  set  of  arcs  in  G  -  G  is  not  empty,  P^t.  includes  successive 

q+1  q  q+J 

sequences  of  arcs  of  ;  let  be  the  last  of  these  sequences  when  P^j+^  *s 

described  from  0  to  D  .  C.  forms  a  route  from  0.  (e  G  )  to  D.  (e  G  )  since 

3  3  q  3  q 

C.  is  a  subset  of  P^,.  . 

j  q+j 

The  part  of  P^ , .  from  0  to  0.  is  included  in  G  for  if  it  were  not,  it 
q+J  J 

would  be  possible,  from  Lemma  2,  to  construct  another  route  shorter  than  tt^  and 
which  would  have  the  arcs  of  C.  outside  of  G  thus  contradicting  Lemma  4  and 

j  q 

similarily  for  to  D  ,  proving  3  a  single  path  . 

Suppose  that  a  proper  subset  of  forms  a  cycle.  If  this  cycle  was  a  0-cycle 

from  Lemma  3,  n  would  be  infinite  and  the  algorithm  would  have  been  terminated  at 
the  end  of  Step  1. 

If  this  cycle  had  a  positive  length,  it  would  be  possible  by  deleting  it  from 
Pj  to  construct  a  route  not  completely  included  in  G^  and  of  length  <  tt^  ,  thus 
contradicting  Lemma  4  again.  Thus  if  there  is  a  cycle,  it  can  only  be  simple. 

q.e.d. 

The  following  lemma  is  not  crucial  for  the  proof  of  the  theorem,  however,  it 
helps  understanding  of  the  algorithm. 


Lemma  6 .  -  tt  >  tt  . 
-  q  q-1 


Proof : 


It  is  impossible  that  tt  <  ti  _  since  P^f.  is  not  among  the  q  shortest 

q  q-1  q+J 


V 

\ 
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routes . 

Suppose  it  =  7i  .  .  Then  3  (x,y)  e  A  Pi  Pq , .  for  which  r  (y)  +  n(y)  = 
q  q-J  q  q+1  q 

ir  =  Tt(j_2  but  then  P^j^  would  have  been  discovered  at  a  previous  passage  and 
(x,y)  would  have  been  added  to  ,  thus  a  contradiction. 

q.e.d. 


Theorem 

The  algorithm  leads  to  a  solution  in  a  finite  number  of  steps,  or  shows  that 
no  solution  exists. 


Proof: 


In  Steps  2  and  3  the  sequence  P^  ,  P^  is  generated  so  that  this 


q+1 


q+k 


sequence  is  at  hand  at  beginning  of  Step  1.  Then  P^+^  is  either  the  shortest 

route  completely  included  in  ,  p|j+^  s  or  the  shortest  route  containing  some  arc 

in  A  i  P^,,  (Lemma  A), 
q  q+1 

If  the  number  of  routes  in  G  ,  (k  +  q)  ,  is  less  then  k  and-  if  there  does 

q  q 

not  exist  routes  outside  those  in  G  (ir  =  °°)  the  problem  is  infeasible. 

q  q 

Each  time  Step  1  is  described,  some  arcs  are  added  to  G^  so  that  we  can  only 
go  a  finite  number  of  times  through  Step  1. 

We  cannot  cycle  indefinitely  in  Step  3  since  each  time  we  go  through  Step  3, 
at  least,  a  new  route  is  added  to  the  updated  sequence  and  we  do  not  want  more  than 
p  routes. 

q 

Now  we  can  go  only  a  finite  number  of  times  through  Step  2  since  at  each 
passage  the  index  j  is  increased  by  one  unit  and  when  j  >  n  we  go  to  Step  1. 

q.e.d. 


A.  Variant  Of  The  Preceding  Algorithm  To  Find  The  k  Shortest  Chains 

The  basic  idea  of  this  variant  is  to  perform  the  algorithm  of  Section  2  and 
to  keep  the  routes  thus  generated  only  if  they  are  chains.  Two  main  differences 


11 


with  the  initial  algorithm  justify  a  rewriting  of  the  algorithm: 

-  When  a  0-cycle  is  met,  no  new  chain  will  ever  be  constructed. 

-  Routes  of  order  u  >_  1  will  not  give  rise  to  any  chain  and  Step  3  can  be 
skipped. 


Step  0* .  -  Find  all  the  shortest  chains  from  0  to  D  :  P|  P*  . 
If  q  ^  k  ,  terminate. 

Determine  for  all  x  e  N  : 

7i (x)  =  shortest  distance  from  0  to  x  . 
ri(x)  =  shortest  distance  from  x  to  D  . 
k'  =  0  . 

q 


Step  l1 .  -  Define: 


(1) 

(2) 


C  (x) 
q 


Min(y,X)  c  Aq‘"(y)  +  a(y'x)) 
7l(x) 


for  x  e  N 

q 


for  x  i  N 

q 


E  *=  {x  e  N  |  e  (x)  +  n(x)  =  Min  [c  (y)+n(y)l  "  ti  } 
q  ycN  **  ^ 


If  ti  =  00  and  k'  +  q  <  k 

q  q 

Determine  m  by: 


terminate,  the  problem  is  infeasible. 


£(P 


«q 

q+m 


) 


<  71 

-  q 


<  £  (P 


»q 

q+m+1 


) 


Let: 


G'  ,  *  G* 

q+m  q 

P » q+m  =  P  t  q 

q+i.rfi  q+:,>i 


q  =  q  +  m 


i  1 . . 


-  m) 


3  =  0 
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For  each  x  e  E  generate  the  routes  R.  .  from  0  to  x  , 

(ox) 

the  last  arc  of  which  belongs  to  and  the  length  of  which  is 

C^(x)  ;  this  is  easily  done  by  backtracking  in  the  definition  of 

C  (x)  and  tt(x)  .  Similarly  generate  all  the  chains  C,  from 
q  (xD) 

x  to  D  of  length  n(x)  . 


By  combining  a  route  R,  N  which  a  chain  C,  in  such  a 

(ox)  (xD) 

way  that  no  newly  constructed  route  from  0  to  D  contains  a 
cycle  counted  more  than  once,  a  finite  number  n  of  routes  from 
0  to  D  of  length  rr^  ,  ^q+1  ^q+n  are  8enerate<^*  Among 

those  are  all  the  chains  of  length  and  having  a  nonempty  inter¬ 

section  with  A 

q 

Go  to  Step  2 ' , 


Step  2 ' . 


j  "  j  +  1 

If  j  >  n  go  to  Step  1' . 
If  Pq,  is  a  chain,  let 

q+j 

Let 


(q  +  1  if  ^q+j  *S  3  c^a^n 

q  otherwise 


If  the  set  of  arcs  in 


G’  U  Pq,, 
q  q+j 


is  empty,  let 


p.q'  =  p.q 

q’+i  q+i 

q  =  q' 


1  <  i  <  k’ 

-  -  q 


Go  to  Step  2*. 
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If  the  set  of  arcs  in  G’ ,  -  G’  is  not  empty,  it  forms  a  simple 

q  q 

route  C.  from  0.  to  D.  .  Let  C!  be  the  subset  of  C.  which 

3  j  1  J  j 

is  a  chain  from  0.  to  D.  . 

J  j 

Let 


R'  =  (R!  ..  R')  ,  T*  =  (T*  ..  T') 
l  r  s 

be  respectively  a  set  of  chains  from  0  to  0,  ,  from  to  D 

obtained  by  taking  parts  of: 


-  Those 

chains 

from 

0 

to  D 

which 

contain 

0.  . 
3 

-  Those 

chains 

from 

0 

to  D 

which 

contain 

D.  . 

’  •  • 
J 


Generate  new  routes  in  G^f  ,  according  to  increasing  lengths,  by 
routes  taking  union  of  a  chain  in  R'  ,  chain  c!  and  a  chain  in 
T*  .  Keep  only  those  routes  which  are  chains  and  insert  them  in 
the  sequence  so  that  the  new  sequence  be  still  ordered 

according  to  increasing  lengths.  Generate  new  chains  until 
exhaustion  of  possibilities  or  until  the  k  -  q  +  j  -  1  first 
elements  of  the  new  sequence  have  been  generated. 

Let 


q  =  q» 

Denote  by  the  elements  of  the  sequence  just  generated 

and  by  k'  their  number. 

q 

Go  to  Step  2 1 . 

This  algorithm,  as  well  as  algorithm  of  Section  2,  can  be 
slightly  modified  in  order  to  produce,  respectively  all  the  chains 
(routes)  which  are  shorter  than  a  given  length  L  (or  within  x% 
of  the  length  of  the  shortest  chain  (route)). 
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5.  Examples 

Suppose  we  want  to  find  the  four  shortest  routes  and  the  four  shortest  chains 
in  the  following  graph: 


a)  Four  Shortest  Routes 

1)  Step  0.  -  The  shortest  route  is  unique  =  [0,  1,  2,  D]  =  {(01),  (12), 

(2D)}  A(PX)  =  6 

2)  Step  1.  - 


Nodes 

71  (x) 

n(x) 

Cj_(x) 

n(x)+c1(x)  , 

0 

0 

6 

3 

9 

1 

1 

5 

00 

OO 

2 

4 

2 

00 

OO 

3 

6 

2 

6 

8 

D 

6 

0 

8 

8 

P^  does  not  exist 
P'  =  [0,  1,  2,  3,  D] 


n  =  1 
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6)  Step  2,  - 


j  =  1 


S  '  G2  U  *3 


i.e. 


4(P l)  -  11 

If  we  go  to  Step  3  and  generate  the  shortest  route  of  order  1»  we 
find 


«=  [ 0-1-2-u- 1-2-0- 1-2 -D] 
i(P 4)  -  12 

3 

"  and  we  insert  it  in  the  sequence  after  . 

Go  to  Step  2. 

7)  Step  2.  -  j  *  2  >  n  -  1  ,  go  to  Step  1. 


8)  Step  1.  - 
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Nodes 

n(x) 

C3(x) 

n(x)+c3(x) 

0 

6 

00 

00 

1 

5 

00 

00 

2 

2 

CO 

0O 

3 

2 

8 

10 

D 

0 

00 

00 

E  =  {3} 

*3  =  11 


and 

P4  =  p3  =  JO- 3-D] 


terminate. 


b)  Four  Shortest  Chains 

Up  to  stage  (6)  the  two  algorithms  give  the  same  result  i.e. 


P{  -  [0,  1,  2,  D] 

■=  (0,  1,  2,  3,  D] 


6')  Step  2 * .  - 


j  -  1 

q'  =  q  =  2 

G2  -  G2  U  53 


No  new  chain  in 


does  not  exist. 


k* 

k2 


=  0  . 


Go  to  Step  2 ' . 


7')  Step  2  * .  -  j  =  2>n  =  l  go  to  Step  1' 


Nodes 


S2(x) 


n(x)+c2(x) 


does  not  exist  so  m 

P3  =  [0,  3,  D] 
n  =  1 


Go  to  Step  2 ’ . 


-  j  -  1 


is  a  chain:  P^ 


-  [0,  3,  D] 


10’) 


G3  "  G’  U  P3  =  G 


There  is  no  chain  outside  P|  P^ 


does  hot  exist,  =  0  . 
Go  to  Step  1' . 


Nodes 


n(x)  I  c0(x)  I  n(x)+c„(x) 


.19 


+  3  =  0  +  3  <4 

Problem  is  infeasible. 


V 

1 
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